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1. Introduction 

The science of graphs and networks is nowadays a well studied topic receiving more 
and more interest due to its large number of applications in the understanding of 
complex systems in, technological, biological and social structures. We refer the 
interested reader to the literature El G] . As a basic starting point in the analysis of 
the cooperative behavior typical of such complex systems, the study of the simplest 
dichotomy-interaction models, such as the Ising model, plays a crucial role and in 
recent years successful efforts have been made in solving this model E] , which, on 
networks, and especially on the so called "scale free networks", may manifest critical 
phenomena very different from those over regular lattice 1 systems G 7 . The 1 Ising 
model with only ferromagnetic couplings, however, allows only a "friendly" interaction 
between the elements of the network, Jb = +J, so that only a possible two-state 
"friendly" collective behavior is allowed. For taking into account that the interaction 
between the elements of the network can also be unfriendly, J5 = ± J, the Ising spin 
glass model should instead be considered [S]. 

In the part I of this work jU], we have introduced a general mapping between a 
given Ising spin glass model and a related Ising one; the related Ising model, which, 
roughly speaking, is the model in which the couplings are uniform. The mapping 
turns out to be of immediate application when the solution of this related Ising model 
is known. 

The mapping becomes exact when the dimensionality of the system goes to 
infinite. More precisely, we have seen that: for a system having a number of first 
neighbors per vertex proportional to D, the mapping applies exactly as D — > 00 
(infinite dimensional in a strict sense) in all the paramagnetic region; whereas, more 
in general, the mapping becomes an exact method to find the upper critical surface 
when, in the thermodynamic limit, the total number of infinitely long paths per vertex 
goes to infinite and, choosing randomly two of them, the probability that they overlap 
each other for an infinite number of bonds goes to zero (infinite dimensional in a large 
sense). 

In |H] we have applied the mapping to several models in which the above condition 
is satisfied, including, generalized Sherrington-Kirkpatrick models, which are infinite 
dimensional in a strict sense, and the spin glass model defined over the Bethe lattice, 
which is infinite dimensional in the large sense. 

In this paper we will apply the mapping to general graphs and networks which 
arc infinite dimensional in the large sense. First, we will consider generalized static 
tree-like structures allowing also for a finite number of loops per vertex. Second, we 
will analyze the case of random graphs where an ensemble of graphs G, together with a 
probability P(g) 1 g G G, are given and all the graphs g are infinite dimensional in the 
large sense, which, in particular, includes the so called "equilibrium random graphs" 
0. In both the cases, we find very general formulae for the exact upper critical surface: 
parmagnetic-spin glass and paramagnetic-disordered ferromagnetic/antiferromagnetic 
transitions. Whereas for the generalized tree-like structures these critical surfaces 
are here for the first time given, we observe that, for the random graphs, the critical 
surfaces have been just derived, but confined to a particular subclass of the equilibrium 
random graphs ;8;. Furthermore, beyond the equilibrium random graphs, we provide 
a formula which applies exactly in the most general case for networks whose bonds 
are each other independent random variables. This formula gives us the upper critical 
surface of the Ising spin glass model in terms of the critical temperature of the Ising 
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model over the network, which, in general, has to be determined numerically, but, as 
a non random model, its measure turns out to be hugely easier than a direct measure 
for the spin glass model. 

The paper is organized as follows. In Sees. 2 and 3, we recall briefly the models 
we are considering and the main result of the Ref. [5]. In Sec. 4, we apply the mapping 
to static generalized tree-like structures. In Sec. 5, a brief introduction for defining 
the random graphs is given. In the subsection 5.1, we recall the critical temperature 
of an Ising model defined over an ensemble of random graphs, which, by using the 
mapping, will be used in the following subsection 5.2 to easily find the upper critical 
surface of a generic Ising spin glass over an ensemble of random graphs. Finally, some 
conclusions are reported in Sec. 6. 



2. Models 



In this paper we will apply the mapping to systems defined over a graph g of N 
vertices whose set of links T g will be defined through the adjacency matrix of the 
graph, g id = 0, 1: 

T g = {b= (i b ,ib) ■ h,jb G .9, gi h ,j h = 1, ib < jb}- (1) 

The fully connected graph will be indicated with Tt: 

T f = {b = (i b ,ib) ■ ib,jb = 1, • ■ ■ ,N, i b < j b }. (2) 

The Hamiltonian of these spin glass models can be written as 

H{{ai};{J b }) = - Jb&b (3) 
ber g 

where the J b 's are quenched couplings, and a b stays for the product of two Ising 
variables, a b = <ii b Oj bl with i b and j b such that b = (i b ,j b ). 
The free energy F is defined by 



~0F^J dV({J b })log(Z({J b })), (4) 
where Z ({ J b }) is the partition function of the quenched system 

Z{{J b }) = Y,e- m{rTih{Jb}} \ (5) 

and dV ({ J b }) is a product measure given in terms of normalized measures dfi b > 
(we are considering a general measure dfx b allowing also for a possible dependence on 
the bonds) 

dV {{J b }) = J] dfXb (Jb) , / d^ b (J b ) = 1. (6) 
ber f 

We take the Boltzmann constant Kb = 1. A generic inverse critical temperature of 
the spin glass model, if any, will be indicated with f3 c - 



3. Mapping 

In the part I of this work we have provided a mapping between an Ising spin glass model 
and a related Ising one. The mapping is exact on and above the critical temperature 
when, in the thermodynamic limit, the dimensionality is infinite. For the aims of the 
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present work, we will need to recall only the following definitions and rules of the 
mapping. 

Given a spin glass model trough Eqs. (11 Kill , we define, on the same set of links 
F s , its related Ising model trough the following Ising Hamiltonian 

ff/(W;W) = -E J ^ ( ? ) 

6er 9 

where the Ising couplings have non random values such that V b 1 b' G T g 

4° = J b ] if d ^b> = d^b, (8) 
4 I] ± if J dfi b (J b )J b ^ or J dpi b {J b )Jl > 0. (9) 

In the following, a suffix / will be referred to the related Ising system with Hamiltonian 
|7|>. Finally, with z b (zjf^) we will indicate the universal parameters z b = tanh(/?Jt>) 
(z^ = tanh(/3J^)), and with w b (wfp) their critical value for (3 C (/?c^). 



3.1. Case of a uniform measure (same disorder for any bond) 

If d/j, b = dfi for any bond b of r g , the related Ising model corresponds to a 
uniform Ising model having a single coupling = and its critical behavior 
will be characterized by, at most, two points urp = tanh(/3^ j( 7 )) > and 

w af = tanh(/3^ J 1 " 1 ^) < 0, if any, where [3y and p\L are the critical ferro- and 
antifcrromagnctic temperatures of the related Ising model, respectively. If r s is infinite 
dimensional in the large sense, the critical inverse temperature of the spin glass model, 
/3 C , is given by 

/? c = min{/?( SG ),/3f/^)}, (10) 
where pi SG ^ and (3^^ AFS> are, respectively, the solutions of the equations, if any 

dfi tanh 2 (/3f G U) =w { p\ (11) 

d^ t^i F ' AF U b ) = w { p /AF: (12) 
where F or AF, in the l.h.s and r.h.s. of Eq. Q12|l. are in correspondence. 



3. 2. Generalization 

In the case of an arbitrary (not uniform) measure dfi b , useful for example for 
anisotropic models, in which we have to consider even a given number of bond 
dependencies, the related Ising model is defined by a set of, typically few, independent 
couplings { J^}, trough Eqs. i18l!)l) and its critical behavior will be fully characterized 
by the solutions of an equation of the type Gi({z b }) = 0. If T g is infinite dimensional 
in the large sense, Eqs. (|1(JI12|) are generalized as follows. The critical inverse 
temperature of the spin glass model (3 C is given by 

/3 c = min{/3p),/3f/^)}, (13) 
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where /3 C SG ' > and /3c A are, respectively, solutions of the two set of equations 

Gj ({ J d^b tanh 2 (/3( 5G) J 6 )}^) =0, (14) 

Gi ({ J d\x h tanh(/3f ' AF ^ J b )]j = 0, (15) 
3.3. Phase diagram 

Equations (|1C)111|I and their generalizations I |13I15|I give the exact critical 
paramagnetic-spin glass (P — SG), Pc , and paramagnetic-F/AF (P — F/AF), 

IT? I A T?\ 

f3c , surfaces, the stability of which depend on which of the two ones is the 
minimum. In the case of a uniform measure, the suffix F and AF stay for ferromagnetic 
and antifcrromagnctic, respectively. In the general case, such a distinction is not 
possible and the symbol F/AF stresses only that the transition is not P — SG. 

4. Spin glass over generalized tree-like graphs (static graphs) 

In the part I of this work [j]], as application of the mapping I|11I12[) . starting from 
the critical value of the universal parameter of the Ising model over a Bethe lattice of 
coordination number q, 

„,M = fftW tV)\ = 1 



= tanh(^J^J = — -, (16) 

we have recovered the critical surface of the corresponding Ising spin glass model: a 

spin glass (SG) and a disordered ferromagnetic (F) transitions take place at /3 C SG '' > 
( 

and (3c , respectively solutions of the two following equations 

dfi tanh 2 (/3f G) Jt) = — -, (17) 

d[i tanh^^Jfc) = — (18) 

which, by using Eq. i|10|) . give the upper phase boundary. 

We now generalize Eqs. I|17|l and H18J) to a spin glass defined over a given 
generalized tree-like graph g. Here g is said to be a generalized tree-like graph if 
for any vertex there is at most a finite number of loops (see Fig. QJ . This application 
provides a new non trivial result. We will exploit the concepts and the conclusions of 
R. Lyons Q3]. 

Let us choose a vertex as the root 0. If v is a vertex, \v\ will indicate the number 
of links of the shortest path connecting to «. A cutset H will indicate a finite set of 
vertices having the property that every infinite path starting from intersects II. In 
other words a cutset generalizes the concept of shell used in the Bethe lattice. Finally, 
to generalize the concept of branching, the branching number of T g , denoted by br r g , 
is introduced: 

br r g = inf J A > : inf A_ ' 1 ' 1 = f • ( 19 ) 
I veil J 

It is easy to check that in the case of a regular lattice of coordination number q, one 
recovers br T g = q — 1. Note however that, in general, br T g does not coincide with 
the average of the number of branches per vertex. 
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In |13| is proved that in the uniform Ising model with positive coupling jV' 
defined over a generalized tree like-graph, the ferromagnetic phase transition occurs 
at 

y W = tanh(^)j«)=^-L-, (20) 

which generalizes Eq. (| 1 f>|) . Correspondingly, by observing that a generalized tree-like 
graph is infinite dimensional in the large sense (see Sec. 3 of [§]), and by using the 
general rule (|11I12|) holding for an Ising spin glass with a uniform measure, dfj, b = d/u,, 
wc immediately find that, a spin glass and a disordered ferromagnetic transitions, take 
place at j3 c and (3 C , respectively solutions of the two following equations 

J dfi tanh 2 (4 SG ) J b ) = ^-L-, (21) 

d M tanh(/?f >J 6 ) = ^r- (22) 

Equations l(2"T|) and l(2"2")) generalize (fTTjl and (fT%|l to tree-like structures. 

We can generalize even more the above result as follows. In ^2] is proved that for 
an arbitrary, i.e. non uniform, Ising model with positive couplings {jj, 1 ^} defined over 
a generalized tree like-graph, the ferromagnetic phase transition occurs at a critical 
temperature given by 

T« = inf { T : inf V TT tanh( J Ht) /T) = \ , (23) 



inf J T : inf £ [] tanh( J b(r) /T) = I 
[ ven t<v J 



where, for any vertex v £ II, the product in the r.h.s. is extended to all vertices 
r belonging to the shortest path from to v (r < v), and b(r) denotes the link 
connecting r and its previous point along the shortest path for v. Formally, in terms 
of the universal quantities {z^} = {tanh(/3J b )}, Eq. (|23|l says that the critical point 
{w^} of this non uniform Ising model is located at 

{wi T) } = sup { {z b } : inf ^ 17 *«r) = \. (24) 



sup I {z b } : inf II Zf >H = f 
I veil t<v I 



We can now, correspondingly, apply the rules I|14ll5f) to a spin glass with a 
measure, dV = Y\ beTf d^b, and we find that a spin glass and a disordered ferromagnetic 

transitions take place at /3c and (i c , respectively solutions of the two following 
equations 

/3p G > = sup J p : inf £ ]J f d^ b{r) tanh 2 (/3J b(T) ) = I , (25) 



I3 ( C F) = sup J p : inf £ J] / d Mb(r) tanh(/?J 6 (r)) = L (26) 
[ uen i-<t> J 

Finally, the upper critical surface follows by applying Eq. (jlOJl . whereas the 
multicritical point P — F — SG is the simultaneous solution of Eqs. 125(1 and (12 Gil 
(system of equations). 
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Figure 1. An example of a generalized tree-like graph with several loops. 

Note that, in the Ref. ^2Jj Eqs. 'I--*' 1 an d (|26() . for a disordered ferromagnetic 
Ising model, were already found but limited to measures dfib with positive support of 
the couplings, { > 0}, whereas Eqs. (|2"T|) and (|2"5|) for the spin glass boundary, were 
still unpublished. 

5. Random graphs 

In the case of random graphs, an ensemble of graphs G, together with a probability 
distribution P(g) for the elements g of G, are given. The free energy of a spin glass 
model defined over the ensemble G is then evaluated as 




where Z g ({Jb}) is the partition function of the quenched system in the graph g with 
couplings {Jb}- 

In the following, we will apply explicitly the mapping to the so called "equilibrium 
random graphs" which arc graphs maximally random under the constraint that the 
degree distribution of a vertex is a given one, p{k); k being the number of connections 
of the vertex. For this family of graphs (but not only) the corresponding sets of 
links T g arc infinite dimensional in the large sense. In fact, even if the number of 
infinitely length loops per vertex is infinite (see, g.e., and 22); due to the fact 
that the bonds are randomly distributed with the only constrain that they have to 
generate the given degree distribution p(k), the probability that two of these loops 
overlap each other for an infinite number of bonds is zero. More in general, the 
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elements of any ensemble G are infinite dimensional in the large sense under the very 
weak condition for the distribution P(g) that the bonds are, each other, independent 
random variables. Rather, the fact that the number of finite loops is infinite in certain 
scale free networks, has been recently raised as a possible doubt about the validity of 
the random tree approximation for the solution of the Ising model on the equilibrium 
random graphs that we are about to recall in the following subsection. 



5.1. Ising model over equilibrium random graphs 

Before to consider the Ising spin glass over a random graph, we need to recall the 
general solution for the critical temperature of the Ising model over equilibrium 
random graphs EJ- Let us indicate with (f(k)) = J2^L P(k)f(k) the average 
of the function degree /(fc) with respect to the degree distribution p(k). The critical 
inverse temperature (3 C of the uniform Ising model with a coupling J is given by the 
following relation 

(k) 

tanh(/3 c J) = {k2)p _\ k)p (28) 

It is interesting to see that this result can be approximately derived even by using 
directly our mapping. Let us consider a random graph with N vertices and L bonds 
so that the average degree is (fc) = 2L/N. A given couple of vertices among the total 
of N(N — l)/2 couples, will be connected by one of the L bonds with a probability 
p = 2L/ (N(N — 1)) = (k) p /(N — 1). Therefore, if we neglect the correlations among 
the bonds, an Ising model with coupling J defined over the graph will have an effective 
distribution of couplings described by 

dV({J b }) = J] dfi b (J b ), (29) 
ber f 

where 

dn(J h ) = 5(J b -J)p + S(J b ){l-p). (30) 

According to Eqs. I|7I9|) . the related Ising model of this disorder model, whose disorder 
is given by Eqs. (|29J) and (|30|l . is the Ising model defined over the fully connected 
graph Tf for which a ferromagnetic or antiferromagnetic phase transition occurs at 

pi^J^N = 1; J^ 7 ' being a generic coupling, positive or negative, respectively. 
Therefore, if we now use the measure l(30() in Eq. 1)12(1 . and Wpj AF = tanh(/3^ ) j( J )) -> 
1/N for N large, we find that the critical inverse temperature j3 c of the Ising system 
over the random graph, is given by tanh(/3 c J) = 1/ (k) , which equals Eq. (|28|l only 

in the case of the classical random graph, for which (fc 2 ) p — (k) p — {k) 2 p [TUj. The 
error is due to having neglected the correlations among the bonds or, in other words, 
in assuming for dV({ J&}) a product measure. On the other hand, in a random graph, 
the distribution of the bonds do not follow the degree distribution of the vertices. In 
fact, as is known, the probability that a bond points to a node of degree k is not p(k), 
but p(k)k/ (k) p . Correspondingly, one expects that, as an effective calculation, in Eq. 
l|5U|l the degree average of a vertex (k) p , should be substituted with the degree average 
of the end of a bond: (k 2 ) p / (k) p . By implementing this substitution one arrives at 

tanh(/3c 7 ' ) J' 7 -*) = (k) / (fc 2 ) , which constitutes a good approximation to the exact 
one G3 HI- 
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5.2. Spin glass over random graphs 

As pointed out above, in a random graph it is not possible to use an effective product 
measure dV({ Jb}) which takes into account rigorously also of the degree distribution. 
However, for an exact treatment we can proceed as follows. Let be given some 
disorder dfi(Jb) (we consider here only the case of a same disorder for any bond). 
We observe that, in Eq. f2"?| . the averages with respect to dV({ Jb}) and to P(g) are 
interchangeable, so that, if for any fixed graph g we apply our mapping, (see Eqs. (67) 
of Ref . jU] ) , for the free energy F we find 

-l3F^N\og(2) + Y J P{9) E / dfx b (J b )log(co S (/3J b )) + Nj2P(9M9), (31) 

where <p(g) is proportional to the non trivial part of the high temperature expansion 
of the free energy density over the graph g and, at least in the limit /3 — > f3~(g), 
it can be evaluated by <p e ff{g) through the mapping equations (24-25) of Ref. jHj, 
j3 c (g) being the critical inverse temperature of the spin glass model defined over the 
(static) graph g and has to be determined according to Eqs. 1)10112(1 . In other words, 
for any fixed graph g, the mapping says that the critical behavior of the system is 
encoded in the effective density <^ e / / (<?) which is nothing else that the non trivial part 
of the high temperature expansion of the uniform Ising model defined over the graph 
g and in which, for any b £ T g , the bond parameter tanh(/3J) is simply substituted 
by J dfi( Jb) tanh(/3J;,) or J d/j,( Jb) tanh 2 (/3Jb), depending on which transition we are 
interested in: disordered ferromagnetic or spin glass, respectively. Since this is true for 
almost any graph g G G (with respect to the measure P(g)), we have that, in the limit 
/? — ► Pcid)' the ensemble G of the Ising spin glass models can be effectively read as an 
ensemble G of uniform Ising models in which the parameters of the high temperature 
expansion are given by J dfi(Jb) tanh(/3J&) or J dn(Jb) tanh 2 (/3Jf,), depending on which 
transition we are interested in. Due to this Ising-representation, by using Eq. (I28|l 
(which holds for an ensemble of uniform Ising models and where now (3 C and J must 
be meant as quantities of the related Ising model, i.e., (3^ and J^\ respectively) in 
Eqs. Ijllfl and i|12|) applied to an ensemble of random graphs with weights P(g) and 
bond parameters J dfi(Jb) tanh(/3Jb) or J dfi(Jb)tanh 2 (pjb), we find immediatly the 
following disordered ferromagnetic or spin glass critical temperatures, respectively 

Ik) 

dv(Jb) tanh(& J b ) = X lp , (32) 
\ K lp ~ \ K )p 

Ik) 

dn{J b ) tanh 2 (/3 c J b ) = ,, 2 , ■ (33) 

)p - W p 

Finally, by applying Eq. I|10|) . the upper critical surface of the phase diagram follows, 
whereas the multicritical point P — F — SG is the simultaneous solution of Eqs. (|32|) 
and l|3l^) (system of equations). 

Equations H32I33J) . together with their consequences, have been just derived in 
the Ref. [5j in the framework of the replica approach, where an analysis of the order 
parameters has also been done. We stress however that, unlike where the random 
graphs considered were confined to the particular class generated by the so called 
"static model" ^H], our derivation shows that Eqs. IpT^ and (|3*3"|) hold for the most 
general equilibrium random graph. More in general, by using the same above scheme, 
given an ensemble of random graphs G distributed with some P(g), if the set T g is 
infinite dimensional in the large sense for almost any g G G with respect to P(g), 
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as happens when the bonds are independent random variables, the critical disordered 
ferromagnetic and spin glass critical temperatures are the solutions of the two following 
equations, respectively 

fd/i(J h ) tanh(/3 c J b ) = wW (G; P) (34) 

J d/i(J 6 )tanh 2 (/? c J&) = #(G;P), (35) 

where w^ 1 ' (G; P) is the (universal) Ising critical parameter in the ensemble G with 
distribution P of the uniform Ising model having coupling J^p; w^(G;P) = 

tanh(/3c 7 ' 1 (G; P) J^ 7 - 1 ), pi (G;P) being the inverse critical temperature. In the case 
of equilibrium random graphs, Eqs. (|34|l and l|35|) reduce to Eqs. I|32|) and i|33|) . 
respectively. 

As anticipated in the introduction, we stress that, in general, a possible direct 
numerical experiment to measure the Ising critical parameter u)W(G;P), turns out 
to be hugely easier rather than a direct measure for the critical temperature /3 C of the 
spin glass model. 

6. Conclusions 

A general mapping presented in the part I of this work |5j , allows to derive the upper 
critical surface of an Ising spin glass model starting from the critical temperature of 
a related Ising model. When the dimensionality of the set of links T is infinite, the 
mapping is exact. In the present work, we have applied such a mapping to the case 
of static generalized tree-like structures and networks. 

Since in a generalized tree-like structure the number of loops per vertex is finite, 
these structures turns out to be infinite dimensional. Therefore, by using the known 
solution for the Ising model over a generalized tree-like structure, we have derived 
the upper critical surface. We stress that, the result found for the ferromagnetic 
disordered transition, holds in the most general case, regardless of the measure d/j,b of 
the disorder, whereas the formulae for the spin glass transition were still unpublished. 

Concerning networks, which are ensembles of random graphs characterized by 
some distribution P(g), the dimensionality is infinite as far as P(g) does not generate 
infinite long paths with an infinite overlapping of bonds, as happens in the so called 
equilibrium random graphs and, more in general, when P{g) is such that the bonds 
turns out to be each other independent random variables. For the case of equilibrium 
random graphs, whose P{g) can be completely described in terms of the only degree 
distribution of the vertices p(k) , the known solution for the Ising model allowed us to 
find exactly the upper critical surface of the Ising spin glass model over the network. 
In this case, as in the Ising model, the critical surface is expressed in terms of the first 
and second moment of p(k). We stress that this result is very general, holding for any 
kind of equilibrium random graph and not only for the so called "static graph" |Hj. Yet, 
we have provided a more general formula which, in terms of the critical temperature 
of the Ising model, applies exactly to any network as far as almost (with respect to 
the measure P(g)) any graph g is infinite dimensional (as happens when the bonds are 
each other independent random variables). Unlike the case of the equilibrium random 
graphs, however, now the critical temperature of the Ising model, and, therefore, of the 
Ising spin glass model as well, will be expressed in terms of more complicated averages 
than the simple moments of the degree distribution p(k). We point out however, that, 
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even if in these cases the analytical knowledge of the critical temperature of the Ising 
model is absent, a direct numerical measure for it, turns out to be hugely easier than 
a direct measure for the critical temperature of the spin glass model. 

By looking back at: Eqs. I|17f) an d GHt for the Bethe lattice case; Eqs. (12111 
and (|2"2"|) or and for the generalized tree-like graphs; Eqs. J2U and 
for the equilibrium random graphs; and Eqs. I|34(l and l|35(l for more general random 
graphs, we recognize a sort of common structure for the critical surface and the critical 
behavior. On the other hand, this is not surprising considering that the related Ising 
models of these models are essentially reduced to an Ising model over some tree-like 
graph, if necessary, distributed according to some distribution whose fluctuations has 
to be taken into account. In this sense, we can say that the mapping establishes a 
sort of universality among different Ising spin glass models. 

We conclude by observing that, as already mentioned in the part I, our mapping 
does not seem to be peculiar of the Ising spin glass models. For more general models, 
a mapping at high temperature "random system" — > "non random system" seems 
possible; so that a known solution of a given non random model can be used to find 
the upper critical surface of the corresponding random model over infinite dimensional 
graphs or infinite dimensional networks. 

Acknowledgments 

This work was supported by DYSONET under NEST/Pathfinder initiative FP6, and 
by the FCT (Portugal) grant SFRH/BPD/24214/2005. I am grateful to A. V. Goltsev 
for useful discussions and to F. Mukhamedov for bringing Ref. to my attention. 

References 

[1] Albert R and Barbasi A L, 2002 Rev. Mod. Phys. 74 47 

[2] Dorogovtsev S N, Mendes J F F, 2003 Evolution of Networks (University Press: Oxford) 
[3] Pastor-Satorras R and Vcspignani A, 2004 Evolution and Structure of the Internet (University 
Press: Cambridge) 

[4] Dorogovtsev S N, Goltsev A V, and Mendes J F F, 2002 Phys. Rev. E 66 016104-1 
[5] Leone M, Vazquez A, Vespignani A, Zecchina R, 2002 Eur. Phys. J. B 28 191-197 
[6] Goltsev A V, Dorogovtsev S N, and Mendes J F F, 2003 Phys. Rev. E 67 026123-1 
[7] Dorogovtsev S N, Goltsev A V, and Mendes J F F, 2005 Phys. Rev. E 72 066130 
[8] Kim D -H, Rodgers G J, Kahng B and Kim D, 2005 Phys. Rev. E 71 056115 
[9] M. Ostilli, "Ising spin glass models versus Ising models: an effective mapping at high temperature 
I. General result", cond- mat/0607498 

[10] Erdos P Rcnyi A, 1959 Publ. Math. Debrecen 6 290 

[11] Bianconi G and Marsili M, 2005 J. Stat. Mech. P06005 

[12] Marinari E and Semerjian G, 2006 J. Stat. Mech. P06019 

[13] Lyons R, 1989 Comm. Math. Phys. 125 337-353 

[14] Such an intuitive argument has been just used in the Refs. [I] and [5], though within different 
procedures. 

[15] Goh K -I, Kahng B, and Kim D, 2001 Phys. Rev. Lett. 87 278701 



